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Abstract 

We consider a new approach to the nucleon-nucleon scattering problem in the framework of 
the higher-derivative formulation of baryon chiral perturbation theory. Starting with a Lorentz- 
invariant form of the effective Lagrangian we work out a new symmetry-preserving framework 
where the leading-order amplitude is calculated by solving renormalizable equations and correc- 
tions are taken into account perturbatively. Analogously to the KSW approach, the (leading) 
renormalization scale dependence to any finite order is absorbed in the redefinition of a finite 
number of parameters of the effective potential at given order. On the other hand, analogously 
to Weinberg's power counting, the one-pion-exchange potential is of leading order and is treated 
non-perturbatively. 

PACS numbers: 11.10.Gh,12.39.Fe,13.75.Cs 
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I. INTRODUCTION 



The few-nucleon sector of baryon chiral perturbation theory was first considered in the 
seminal papers by Weinberg [ll H] which triggered an intensive research activity starting with 
the pioneering work of Ref. [3|. In this approach for processes involving > 1 nucleons 
the power counting is applied to the effective potential, which is defined as the sum of 
all A^-nucleon-irreducible diagrams. The scattering amplitudes are obtained by solving the 
Lippmann-Schwinger (LS) or Schrodinger equation. For recent reviews and references see 
e.g. Refs. 

The implementation of Weinberg's program encountered problems. The NN potential 
of the effective field theory (EFT) is non-renormalizable in the traditional sense already at 
leading order. The iteration of the potential in the LS equation generates divergent terms 
with structures which are not included in the original potential. Therefore, the renormal- 
ization of the solution to the LS equation requires the contributions of an infinite number 
of higher-order counter-terms. Here, the infinite number refers to both the loop and the 
chiral expansions. The freedom of choosing the finite parts of these counter-terms is com- 
pensated by the running of the corresponding renormalized couplings. It has been argued 
that the coefficients of the divergent parts of the counter-terms contributing in low-order 
calculations would set the scale of the corresponding renormalized couplings. As a con- 
sequence, even if these couplings were natural at some value of the renormalization scale, 
they would become unnaturally large for slightly different values of this parameter. This 
problem (also analyzed non-perturbatively), in different variations, has been addressed as 
the inconsistency of Weinberg's approach, and an alternative power counting scheme, known 
as the KSW approach, has been suggested [gI, 0, B 0] . In the KSW scheme the troublesome 
one-pion-exchange (OPE) potential is shifted from leading order (LO) to next-to-leading 
order (NLO). The leading-order equation becomes renormalizable, i.e. all divergences can 
be absorbed in the redefinition of the parameters of the LO potential, perturbatively as 
well as non-perturbatively. Moreover, the LO equation is exactly solvable and dimensional 
regularization can be applied. Corrections are treated perturbatively in the KSW approach. 
This guarantees that all divergences are absorbed in parameters contributing to the given 
order. As a consequence, the KSW approach is free of "consistency problems." Unfortu- 
nately this scheme suffers from convergence problems. As the OPE potential is shifted to 
NLO, it is treated perturbatively. It has been argued that the perturbative treatment of 
pions in NN scattering poses problems of convergence [l^, 11, 12]. The results of Ref. 13 



made it explicit that perturbative series do not converge in the KSW approach, at least in 
the triplet channel. 

As a next step, in Ref. [3] a non-perturbative investigation of the renormalization prob- 
lem has been carried out and an expansion of nuclear forces about the chiral limit has been 
suggested. It has been argued that this expansion is formally consistent and is equivalent 
to the KSW power counting in the ^5*0 channel and to Weinberg's power count iiig in the 
^5*1 — ^ Di coupled channels. Higher partial waves have been investigated in Ref. [l5| with 
the result that the perturbative treatment of the OPE potential is not sufficient for a finite 
number of partial waves, e.g. ^Pq and F2. On the other hand, the iterated OPE 

potential produces a cutoff dependence in all waves where the tensor force is attractive. 
This requires additional counter-terms in each of these partial waves, clearly upsetting the 
power counting. The authors of Ref. [l5| conjecture that the "mixture of perturbative treat- 
ment of higher partial waves, resummation of lower partial waves, and promotion of a finite 
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number of counter-terms is the most consistent approach" to the chiral perturbation theory 
for nuclear systems while sub-leading interactions are treated perturbatively. 



There exists an alternative point of view that the " consistency problem" of Weinberg's 
approach is irrelevant and one can perform the calculations in cutoff effective field theory 
(EFT) by suitably choosing the cutoff parameter. According to this point of view, taking 
large values of the cutoff parameter without including the contributions of all (perturba- 
tively) relevant counter-terms is not a legitimate procedure. Instead one has to choose the 
cutoff parameter of the order of the relevant large scale, presumably the mass of the p meson 

0, m H H H in, [2I, 



In the present work we develop a new approach which combines aspects of both points 
of view outlined above. In particular, in our new framework the leading-order amplitude 
is calculated by solving integral equations and corrections are taken into account perturba- 
tively. Analogously to the KSW approach, the (leading) renormalization scale dependence 
to any given finite order is absorbed in the redefinition of a finite number of parameters of 
the effective potential at the given order. On the other hand, analogously to Weinberg's 
scheme, the OPE potential is of leading order and is treated non-perturbatively. 

In the two-nucleon sector the path integral for the generating functional of the Green's 
functions is not Gaussian and therefore the standard heavy-baryon reduction can not be 
carried out. While the non-relativistic effective Lagrangian can be fixed by matching to 
the (explicitly) Lorentz-invariant Lagrangian, we prefer to work in terms of the latter. We 
start with the higher- derivative formulation of the manifestly Lorentz-invariant Lagrangian 
of Ref. [23] and perform the non-relativistic reduction at the level of Feynman diagrams. 

Our work is organized as follows: In Sec. [TT] the effective Lagrangian in its manifestly 
Lorentz-invariant form and field redefinitions leading to the higher-derivative formulation 
are considered. Section IIIII addresses the issues of renormalization and power counting. In 
Sec. |IV]the NN scattering equations are derived and analyzed. A summary is given in Sec. El 



II. THE EFFECTIVE LAGRANGIAN 



In this section we consider the effective Lagrangian needed for the purposes of this work. 
The standard effective Lagrangian is organized in a derivative and quark-mass expansion 
and consists of the sum of the purely mesonic, the vrN, NN, . . . parts. 



cflf 



The lowest-order mesonic Lagrangian reads [25 



-Tr 



D^U (D^U)^ + — Tr {xU^ + Ux^) , 



(1) 



(2) 



where U is a unimodular unitary (2 x 2) matrix containing the pion fields. The covariant 
derivative is defined as 



and 



••Hi 



X = 2B{s + ip). 



3 



Here, v^, a^, s, and p are external vector, axial-vector, scalar, and pseudo-scalar sources, 
respectively. In Eq. ([2]), F denotes the pion-decay constant in the chiral limit: F.„- = F[l + 
0{rh)] = 92.4 MeV. We consider the isospin-symmetric limit = = m, and the lowest- 
order expression for the squared pion mass is = 2Bm, where B is related to the quark 



condensate {qq)o in the chiral limit [25 



The nucleon field is described with two four-component Dirac fields, 

n 



and the lowest-order Lagrangian of the one-nucleon sector reads [26| 

= ^ {^1,D^' -m+^gA l,l,u^^ (3) 

where = {d^ -|- — iv^i^^)'^ denotes the covariant derivative with the iso-scalar vector 
source v\^\ and the quantities m, m^, F^, and x± (and their derivatives) are defined as 

= U, u^ = iu^D^Uu\ = ^ [u^d^u + ud^u^ - i{u^r^u + ul^u')] , 

X± = u^XU^ ± ux^u. 

o 

In Eq. ([3]), m and stand for the chiral limit of the physical nucleon mass and the axial- 
vector coupling constant, respectively. 

The most general chirally-invariant leading-order NN Lagrangian reads (for the power 
counting see below) 
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a=0 



+ C\y ^r'^757^$ ^r'^757'^$ + ^rV^^^ ^rV^'^^ 



(4) 



where r° is the unit (isospin) matrix and r* [i = 1, 2, 3) are (isospin) Pauli matrices and 
Cj = C] = Cf for all /. Note that the structure '^'75^' is assigned leading order. 

The most general effective Lagrangian includes all terms which are consistent with the 
symmetries of the underlying fundamental theory, i.e. QCD. The above given "canonical" 
form of the effective Lagrangian is obtained from the most general form by using suitable 
field redefinitions jl^]. This is done for convenience. Any two forms of the effective La- 
grangian obtained from each other by using field redefinitions are completely equivalent 
when physical quantities are considered. Note that the coupling constants corresponding 
to the same operator structures in two forms of the effective Lagrangian are in general not 
equal [2^. Switching from one formulation to the other, one needs to take into account the 
new contributions (in general in all operator structures) generated by field transformations. 
There is an infinite number of possible field transformations and correspondingly an infinite 
number of forms of the effective Lagrangian. While all these forms lead to the same expres- 
sions for physical quantities (with certain care), in practical (perturbative) calculations one 
form of the effective Lagrangian could be considerably more convenient than the other. 
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Starting with the above given " canonical" effective Lagrangian one encounters ultraviolet 
divergences in the equations for the few-nucleon scattering amplitudes. This problem can 
be handled by considering a nucleon propagator with improved ultraviolet behavior ^]. To 
that end we employ a "non-canonical" form of the effective Lagrangian generated by using 
suitable field transformations. To obtain the modified nucleon propagator we perform the 
following field transformations 

2N 



^ ^ ^<l> (^L^ -D j , (5) 
where ^ D = ( D ^ 1 7 and obtain 

£S' = ^ $ -^d') (^z7^D^ -m + ^QA TmTs^'^) $ (^', -^') * ■ (6) 

In Eq. ([6]) L is a free parameter, D stands for the spacial components of the covariant 
derivative and $ is some appropriately chosen real function. Possible choices of $ are, e.g., 

$ ^) = (1 + ;^ ^) or $ (L^, X) = exp (-1 . (7) 

The Feynman propagator corresponding to the Lagrangian of Eq. reads 

^Sf{p) = ^ '-—-FiL^p'), (8) 

p — m + te 



with 
respectively. 

Depending on the order of the performed calculations we choose the modified propagator 
such that all loop integrations over the momenta of the intermediate purely nucleonic states 
are rendered finite. We apply dimensional regularization to the resulting effective theory, 
expand the regularized diagrams in powers of n — 4 and subtract 1/ (n — 4) pole-terms, where 
n is the space-time dimension. There is a finite number of pole-terms to any given (finite) 
order in the chiral expansion of physical quantities in the vacuum and one-nucleon sectors, 
and the effective potentials in the two-nucleon sector. No further divergences occur (for 
finite parameters L) neither in the vacuum and one-nucleon sector nor in the equations of 
the two-nucleon sector. Therefore, we can take n = 4 in the equations of the two-nucleon 
sector after renormalizing the effective potential and the corrections to the two-nucleon 
propagator. Note that the field redefinition of Eq. introduces additional interaction 
terms in the effective Lagrangian. Their contributions are important for preserving the 
symmetries of the theory. In practice it is convenient to choose the exponential form of the 
field redefinition with sufficiently large Ni^, [s'j where the additional interaction terms start 
contributing at orders higher than the order of accuracy of the performed calculations. 

The complete dependence of the physical quantities on the parameter L can be absorbed 
in the redefinition of (an infinite number of) coupling constants. Therefore different choices of 
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the numerical value of this parameter correspond to different renormalization schemes. Note 
that the applied regularization scheme is the dimensional regularization. The parameter L 
should not be confused with the cutoff parameter and it does not have to be taken to infinity. 



III. POWER COUNTING 



It is understood that the power counting specified below is satisfied by renormalized 
diagrams, provided the renormalization points are chosen appropriately. In particular, fol- 
lowing Ref. [2I one should choose the renormalization points of the order of small external 
momenta. 



A. Effective potential 



Interaction terms, generated by the one-nucleon sector of the effective Lagrangian 2^ 
are assigned the standard orders of BChPT 30(]. In the two-nucleon sector each derivative 



acting on the pion field as well as each factor of the pion mass is assigned order g^, where 
q denotes a small quantity like the pion mass or a nucleon three momentum. Contracted 
derivatives acting on nucleon fields lead to scalar products of large nucleon momenta, e.g. 
Pi ■ p2- We can rewrite 

2pi ■ P2 = (pj - vr?) + {p\ - m^) - (pi - p2f + 2m^ (10) 

As in the one-nucleon sector, {pj — w?) is counted as g^, while the difference (pi —P2Y is of 
order q^. We arrange the Lagrangian such that terms involving derivatives are accompanied 
by interaction terms where each pair of contracted derivatives is replaced with a factor of 
"+m^" or "— m^". Therefore each pair of contracted derivatives acting on nucleon fields is 
assigned order (we choose the renormalization scheme such that this power counting also 
holds for loop diagrams.) For example, the term ^D^j^^D^^! is rewritten as 

[^D^^^D^"-^ + ^^^^^^) - ^^^^ , (11) 

where the first and second terms are part of the NLO and LO NN Lagrangian, respectively. 

We assign to each diagram contributing to the effective potential a chiral order D de- 
termined as follows. Vertices from the /c-th order Lagrangian count as order q^ , a loop 
integration in n dimensions as g", each pion propagator as and each nucleon propagator 
as g~^, respectively. Note that the terms of the fc-th order NN Lagrangian can only generate 
contributions to effective potential of order k and higher. To any finite order there is only 
a finite number of terms in the two-nucleon sector of the effective Lagrangian. Correspond- 
ingly the effective potential contains a finite number of terms to any finite order. We drop 
nucleon loops, including their contributions in the redefinition of the coupling constants of 
the effective Lagrangian. 

B. Two-nucleon propagator 

We rewrite the fermion propagator as 

Sf{p) = ^——F{L ,p 

p^ — 171'' + ze 
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4'^ (p) + 4'^ (p) + 4^ (p) + 4^ (p) , (12) 



where the projector P+ is defined as 



/•. = ii^. (13) 



with V = (1,0,0,0). The poles in Sp and Sp correspond to particle states, and those in 
Sp^ and Sp^ to anti-particles, respectively. We assign order to Sp\ order g° to the 
second and third terms Sp^ and Sp\ and order to Sp\ 

The undressed two-nucleon propagator is defined as the product of two undressed one- 
nucleon propagators. Inserting the expansion of Eq. (fT2l) we obtain a sum of terms of 
different orders. For nucleon-nucleon scattering the products of Sp^ and Sp^ are enhanced 
in purely two-nucleon intermediate states [H, [Ij. Therefore, we assign the following orders 
to the various terms in the undressed two-nucleon propagator: 

4)4)~g-, 4)4) 



The first term in Eq. (1141) defines the leading-order two-nucleon propagator Go, 

Go = -^SPs?. (15) 

The last four terms in Eq. (|T4l) possess poles on the same side of the real axis in the com- 
plex plane of the zeroth component of the integration momenta. Therefore, if the effective 
potential does not contain singularities in the complex energy plane, they give vanishing 
contributions. 

IV. NN SCATTERING EQUATION 

Defining the two-nucleon-irreducible diagrams as the effective potential V the equation 
for the amplitude (amputated Green's function) reads 

Txa,^,u{P3,P'l,Pl,P2) = Vxa,fiu{P3,P4,Pl,P2) 

/d^k 

xSF,af3iP3 - k) Sf,^5{P4 + k) Tps^^^ (pa -k,p4 + k,pi,p2) , (16) 
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where the Greek indices labehng the components of the Dirac fermion field run from 1 to 4. 
The corresponding physical amplitude is obtained as 

\l Z^pips) \J Z^{pi) ux{pz)ua{Pi) Txa,^u {Ps, P-i, Pi, P2) (pi (^2) \/ Z^{pi) \J Z^{p2) , (17) 

where is the residue of the nucleon propagator. The residue corresponding to the prop- 
agator of Eqs. ([H]) and ffT^ is 

ZM=F{L\f). (18) 



For convenience in the following derivations we write Eq. (1161) symbolically as 
or 

T = V + VGT . (20) 
Next, we expand T, V , and G in small parameters, 

G = Go + Gi + + ■ ■ ■ , (21) 
V = Vo + Vi + V2 + --- , (22) 
T = To + Ti + T2 + ■ ■ ■ , (23) 

where the quantities Xj (X = T, G, V) are of order C(g*). We then substitute the expansions 
in Eq. ( 120|) . and solve T order by order. At leading order we obtain 

To = yo + K)GoTo. (24) 

To determine the physical amplitude order by order we expand the solutions to the Dirac 
equation in small quantities 

'A — p ■ V 

u = Mo + Ml + M2 H = P+ u(p) H P+ u{p) H , 

2m 

u = uq + u^ + u2 + ■ ■■ = u{p) P+ + u{p) p j' + ■ ■ ■ . (25) 

Substituting Eqs. (123!) and fl25|) in Eq. (fT7|) we obtain the expansion of the physical amplitude 
as 

+UqUq Ti uqUq + UqUq To UiUq + UqUq Tq UqUi H ] ^/Z^ ^/Z^ . (26) 

Taking into account Eq. (125|1 we see that to calculate the leading-order on-mass-shell am- 
phtude we need the quantity To = P4.P+T0P+P4.. Applying the projectors P+ to Eq. (p4|) 
and using (see Eq. (fT5l) ) 

Go = G0P+P+ = P+P+Go = P+P+G0P+P+ (27) 

we obtain 

P+P+ To P+P+ = P+P+ P+P+ + V^o Go P+P+ To P+P+ , (28) 
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which is an integral equation for the projected amplitude Tq, 

To = Vo + VoGoTo, 

where Vq = P+P+ Vq P+P+. 

To calculate Tq from the solution to Eq. fl29|) we rewrite Eq. fl2^ as 

To = K) + K) G K) + K) Go To Go K) , 
which, using Eq. fl27j) . takes the form 

To = \/o + H Go + Go To Go . 



(29) 



(30) 



(31) 



The leading-order amplitude To calculated from Eq. can then be used in calculations 
of the next-to- leading-order amplitude Ti. The equation for Ti 



T^ = Vi + GoTo + VoGiTo + VoGoTi 
can be solved exactly. The solution reads 

Ti = K+ToGoVi + V^iGoTo + ToGol^iGoTo + ToGiTo. 



(32) 



(33) 



Analogously to the above case the leading-order amplitude Tq and the next-to-leading- 
order amplitude Ti can be used to calculate the next-to-next-to-leading-order amplitude T2 
etc. 



A. Leading-order equation 

In the following, we will consider the scattering amplitude in the center-of-mass frame, 



where pi = {^/m?l^,p), p2 = {\/ m?' +p'^, -p), p^ = {sj: 



+ p' '^,p' ), 



m? + p' "^^—p' ). Let P = Pi + P2 = (2po; 0) denote the total four-momentum of the 

scattered nucleons where po = \/p'^ + m?. 

The leading-order effective potential consists of the contact interaction part and the one- 
pion-exchange diagram. The leading-order contact interaction potential reads 



•T/«3i4,n«2 



a=0 



ly^a / a a a a \ 



(757a)A^ (757°)^^ ^fsia^ii " i^5la)xf, 



T- ■ T- ■ 



T- ■ T- ■ 



(34) 
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The projected potential is 

a=0 

b=l ^ 



with l,s,m,n = 1,2. Note that the pseudoscalar piece of the leading-order NN Lagrangian 
starts contributing at higher orders. 

The covariant form of the one-pion-exchange potential reads 



« : — : r,-„,, r, 

a=l 



4F2 (p^_p3)2_M2 '''''''' 



• _9a_ [75(jZ<l-^4)]^^[75 {h-^4)]xu a a (nr.. 



a=l 



In order to simplify some technical aspects of solving the leading-order equation we separate 
the pion denominator as 

1 1/1 1 \ 

+ \ T^^ + ^r-—-] . (37) 



and take the first term on the right-hand side of Eq. ( 137|) as the leading-order contribution 
in the effective potential.^ The second term is of higher order and is taken into account 
perturbatively. The resulting OPE potential reads 



^0,TT;Xa,tiu\P ^P) — / J 



cry 

I -^-^ — : r,- 



. _9a_ [75 ih-^A)l^,b5 (^1-^4)] A. , , , . 

^'4F2 (p-_p-)2 + M2 W^n- l-i^J 



The projection reduces Eq. fl38l) to the following form, 

,-T7i3i4,ni2 ^^ _ V- , ^ [^-iP-Pllim [^-{P-P')] 

'' ^OTT;ls,mn\P ^P) " /I ^2 -a2 , 



a=l 
3 



4F2 (p'_p'f + M2 

[a- (P + P')].m[^- (P + P')]/n « 



• _yA_ -yP^P )\sm -yP^P )\ln a a 
^ 4F2 (p'+p'/)2 + ]Vf2 *3i2*4n- V ) 



The leading-order scattering amplitude satisfies the equation 

^S:^ iP3^p^^pi^P2) = ip'^p) - ^/ (^''-^ 



^ Without this simphfication the Bethe-Salpeter equation has to be solved, which is technically more in- 
volved. 
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(a) (6) 
FIG. 1: Two examples of iterations of the LO equation. 



X 



m2 F^{L^,k^) 



(40) 



where 

K, = V^,,c + Vo,. ■ (41) 

As the leading-order potential does not depend on k^, the integration over this variable can 
be performed easily by iterating Eq. (HOl) . performing the ko integration in each term and 
summing up the resulting terms. Doing so we obtain 



-'-0;ls,mn \P iP) ~ ^Q-ls,mn \P iP) 



m 



d^k 



X 



2 J {2 7t) 
F\L\P) 



3 ^0;ls,ab \P ' 



^2 _|_ j I _ ^2 _|_ ^2 _|_ ^ g 



0;ab,mn 



-k,p 



(42) 



Next we split the two-nucleon propagator into two parts, 

1 2m 



^2 _|_ ^2 



(^Po - VP 



+ + ie) (fc^ + m^) ip'^ — + i e 



+ 



+ VA;^ + j —2m 

^2 _|_ j^2 j I p'2 _ p _|_ ^ g 



(43) 



and take the first term as the leading order propagator while including the second part in 
perturbative corrections. The resulting leading-order equation reads 



d^k 



(2vr) 



yiaU^ij ( ft' Z 



FHL\k^^ 



X 



_|_ ^2 j I ^"2 _ /^2 _|_ ^ g 



Q;ab,mn 



-k,p 



(44) 



Let us investigate the "consistency problem" which manifests itself in the sensitivity of 
the scattering amplitude to the choice of the parameter L. For that purpose we vary L from 
750 MeV to infinity. Divergences which occur in the limit L —>■ oo will be referred to as UV 
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divergences. Simple UV counting shows that iterations of Eq. (jUj) lead to finite diagrams 
for L — > oo. For example, the diagrams in Fig. [T]are finite. This differs from the standard 
non-relativistic formulation, in which Fig. [1] (a) contains an overall logarithmic divergence 
~ M'^ m? and the diagram in Fig. [T] (b) contains a logarithmic divergence ~ {QmY . The 
divergences in these diagrams can be absorbed in contact interaction terms of the second 
and sixth order (and higher-order iterations require higher-order interaction terms, up to 
infinity) 

Equation (jH]) can be rewritten in a more familiar form. To that end we multiply Eq. 



with 



m^F{L^,p'^)F{L^,p^) 



[(m2 +p2)(^2 +p'/2)]l/2 

and obtain 

■^0;ls,mn \P ^P) — '^0-ls,mn \P ^Pj 



where 

%{p',P) 



fo(p',p) F{L^,p'^)F{L\p^)m 
[(m^ + p' 2)(m2 -|- p2^ji/2 

Voip',p) F{L\p'^)F{L\p 



2\ ^2 



Vo(p",p) = r/ 2 ^ ^,2\l 2 I -2^11/2 " i^^) 

Note that the potential Voip, q) has a milder ultraviolet behavior than the standard potential 
^o{p',p) due to the denominator in Eq. (H6!) . In the limit L — oo the factors F{L'^,p'^) 
and F{L'^,p'^) approach 1. The OPE part of the effective potential Vo{p',p) of Eq. (H6|) 
then behaves like ~ 1/p^ for large p in the ^5*0 wave (after absorbing the local part in the 
redefinition of the contact interaction term) and generates no divergeiices (perturbatively 
as well as non-perturbatively). This is in contrast to the OPE part of Vq {p',p) which goes 
as ~ The most singular part of the effective OPE potential in the ^5*1 wave requires 

a single counter-term, which is present in the contact interaction part of our LO potential. 
For higher partial waves the singular behavior of the potential is screened by the angular 
momentum barrier and therefore no counter-terms are required. Analyzing the convergence 
criterium for Neumann series we have checked that the integral equations for higher partial 
waves have unique solutions, i.e. the angular momentum barrier provides a repulsion strong 
enough for our problem. The numerical analysis confirms the above conclusions. Some of 
these phase shifts calculated for different values of the parameter L are shown in Figs. [H [3] 
andm Note that different choices of the parameter L correspond to different renormalization 
schemes. By choosing the renormalization scheme properly one can improve the convergence 
of perturbative series for physical quantities. Our figures suggest that the best description 
of the phase shifts at leading order is obtained for L ~ 750 MeV. This indicates that the 
higher-order corrections are smaller for this choice of the renormalization condition. 

We conclude that all counter-terms required for the absorbtion of the strong L dependence 
in the LO equation (H2l) . perturbatively as well as non-perturbatively, are present in the LO 
potential. Therefore, Eq. (H2l) is free of the "problems of consistency of Weinberg's power 
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FIG. 2: ^Sq partial wave np phase shifts for different values of L parameter compared with the 
data points from the Nijmegen PWA [3l[ |. The curves are shown for L = 750, 1000, 2000, 6000, 
20000, 30000, 100000 (MeV), the sohd line corresponding to the largest value. 
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FIG. 3: ^Si partial wave np phase shifts for different values of L parameter compared with the data 
points from the Nijmegen PWA [3l|. The curves are shown for L = 750, 1000, 2000, 6000,10000, 
20000, 50000 (MeV), the solid line corresponding to the largest value. 



counting" encountered in solving the corresponding Lippmann-Schwinger equation 

Note that Eq. fH21) in comparison with the LS equation contains the re- summation of 
an infinite number of higher-order terms. In the authors' opinion both approaches are 
equivalent and self-consistent, provided that one includes compensating terms in the effective 
Lagrangian and chooses the cutoff parameter appropriately in the LS approach. Therefore 
we expect the numerical results for physical quantities in our approach to be comparable to 
that of the cutoff LS equation approach [5|. 



^ The Lippmann-Schwinger equation corresponds to the leading-order term in the 1 /m expansion of Eq. (|42p . 
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FIG. 4: ^Pq partial wave np phase shifts for different values of L parameter compared with the 
data points from the Nijmegen PWA [3ll ]. The curves are shown for L = 750, 1000, 2000, 6000, 
20000, 30000, 100000, 500000 (MeV), the solid line corresponding to the largest value. 



B. Higher-order Corrections 

Corrections to the leading-order amplitude are calculated perturbatively and to any finite 
order can be expressed in a closed form as a sum of a finite number of terms analogous to 
Eq. fl33l) . These expressions can be interpreted as sums of a finite number of Feynmann 
diagrams where the leading-order amplitude Tq is also interpreted as an NN vertex. The 
straightforward analysis of Eqs. fH^ - fj45l) shows that (for L —>■ oo) To(p, j?) cannot diverge 
stronger than ~ Inp (~ Inp') for p — > oo (p' — > cxd). This UV behavior of Tq guarantees 
that all divergences generated in the limit L ^ oo in Feynman diagrams contributing to 
corrections at any finite order can be canceled by counter-term contributions also present 
at the given order. As a result the new semi-relativistic approach suggested in this work is 
free of " consistency problems" inherent in the standard non-relativistic approach. 



V. SUMMARY 



We consider a new approach to the nucleon-nucleon scattering problem in baryon chiral 
perturbation theory. Starting with a Lorentz-invariant form of the effective Lagrangian with 
higher derivative terms we work out a new framework where the nucleon-nucleon scattering 
amplitude is calculated order-by-order in small parameters like the pion mass and nucleon 
three-momenta. We develop a systematic power counting for the effective potential which is 
defined as the sum of all two-nucleon irreducible diagrams. The leading-order amplitude is 
obtained by solving an equation with the LO potential consisting of the contact interaction 
plus OPE. Higher-order corrections to the potential as well as to the two-nucleon propagator 
are taken into account perturbatively. The parameter L, introduced through terms with 
higher derivatives, parameterizes the renormalization scheme dependence of the scattering 
amplitude present in our calculations. Due to the UV behavior of our semi-relativistic two- 
nucleon propagator, analogously to the KSW approach, the strong L-dependence of physical 
quantities to any finite order is absorbed in the redefinition of a finite number of parameters 
of the effective potential at given order. On the other hand, analogously to Weinberg's 
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power counting, the OPE potential is of leading order and is treated non-perturbatively. 
Our new approach preserves all underlying symmetries to the given order of accuracy and 
the self-consistency of the power counting is explicit in this new formulation. 
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